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ABSTRACT 



A major goal of the hydraulic ram survivability program 
for aircraft fuel tanks is the development of an economical 
computer code which accurately predicts the pressure in the 
fluid, the response of the tank walls, and the crack damage 
to the tank due to a ballistic projectile. This study com- 
pares the tank wall response predicted by the recently modi- 
fied SATANS code with the data from studies performed by 
California Research and Technology, Inc., the Naval Weapons 
Center, and the University of Dayton Research Institute. Good 
agreement among the results is obtained for this phase of the 
study. Additionally, a comparison of the predicted crack 
length in the entry wall of the tank using Fahrenkrog's cri- 
terion and SATANS ' predicted circumferential stresses is made 
with actual crack lengths from a test performed by the Univer- 
sity of Dayton Research Institute. The predicted crack lengths 
were in excellent agreement with the actual crack lengths. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 8 

II. COMPUTER STUDIES 11 

A. BACKGROUND 11 

B. ENTRY WALL RESPONSE-A COMPARISON WITH 

CRT RESULTS 12 

C. EXIT WALL RESPONSE-A COMPARISON WITH 

HRP-V1 MOD RESULTS 16 

D. CRACK DAMAGE STUDY 17 

III. CONCLUSIONS AND RECOMMENDATIONS 20 

A. ENTRY WALL RESPONSE 20 

B. EXIT WALL RESPONSE 20 

C. CRACK DAMAGE STUDY 21 

FIGURES 22 

LIST OF REFERENCES 4 4 

INITIAL DISTRIBUTION LIST 45 



5 



LIST OF FIGURES 



1. Finite Difference Mesh 22 

2. NONSAP Finite Element Model of Entry Panel for 

CRT Study 23 

3. Conditions for CRT Calculations and Sphere 

Ballistic Data 24 

4-a. Spherical Projectile Displacement vs. Time 25 

4- b. Spherical Projectile Velocity vs. Time 25 

5- a. Entry Wall Displacement vs. Distance from Impact 

Point at t=200 ysec 26 

5-b. Entry Wall Displacement vs. Distance from Impact 

Point at t=400 ysec. 27 

5-c. Entry Wall Displacement vs. Distance from Impact 

Point at t=600 ysec. 28 

5-d. Entry Wall Displacement vs. Distance from Impact 

Point at t=800 ysec. 29 

5-e. Entry Wall Displacement vs. Distance from Impact 

Point at t=1000 ysec. 30 

5-f. Entry Wall Displacement vs. Distance from Impact 

Point at t=1200 ysec. 31 

5-g. Entry Wall Displacement vs. Distance from Impact 

Point at t=1400 ysec. 32 

5-h. Entry Wall Displacement vs. Distance from Impact 

Point at t=1600 ysec. 33 

5-i. Entry Wall Displacement vs. Distance from Impact 

Point at t=1800 ysec. 34 

5- j. Entry Wall Displacement vs. Distance from Impact 

Point at t=2000 ysec. 35 

6- a. Circumf erential Stress vs. Radius at t=50 ysec. 36 

6-b. Circumferential Stress vs. Radius at t=100 ysec. 37 

6-c. Circumferential Stress vs. Radius at t=150 ysec. 38 

6-d. Circumferential Stress vs. Radius at t=200 ysec. 39 



6 



7. Exit Wall center Plate Deflection vs. Time 40 

8. SATAIJS and HRP-V1 MOD Plate Deflection vs. 

Radius at t=2000 sec. 41 

9. Circumferential Stress vs. Distance from Impact 

Point for crack growth 4 2 

10. Maximum Circumferential Stress vs. Radius for 

crack growth 4 3 



7 



I. 



INTRODUCTION 



Fuel tanks in combat aircraft are probably the most vul- 
nerable system in the aircraft. When a bullet or missile 
fragment penetrates the wall of a fuel tank and enters into 
the fuel, an intense pressure field propagates throughout the 
fuel and impinges on the walls of the tank. This is called 
the hydraulic, or hydrodynamic, ram phenomenon. It is a com- 
plicated process, and a considerable amount of research effort 
has been devoted to it by many investigators in their efforts 
not only to gain an understanding of the phenomenon but also 
to develop an economical computer program which accurately 
predicts the pressure in the fluid and the response of the 
tank walls. Many of the previous studies were presented at 
the Hydrodynamic Ram Seminar held at the University of Dayton, 
Ohio, October 1976. A summary of the hydraulic ram research 
performed at the Naval Postgraduate School is contained in the 
symposium proceedings. Ref. 1, and also in Refs. 2 and 3. 

The computer studies of hydraulic ram conducted at the 
Naval Postgraduate School have utilized a computer code called 
SATANS. A description of the code and subroutines is contained 
in the user's manual. Ref. 4. SATANS is a structural response 
digital computer code which has been used to predict the dis- 
placements, stresses, and strains in a circular wall of a 
fluid-containing tank as a function of time. The fluid-struc- 
ture interaction at the tank wall was accounted for in SATANS 
by the plane wave piston theory. The piston theory was also 
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used in the development of the finite element code BR-1HR, 
which can also be used in fluid-structure interaction prob- 
lems. However, recent studies have shown that the piston 
theory assumption that the transverse pressure variations in 
the fluid can be neglected leads to extremely inaccurate wall 
response predictions except at very early times. Ref. 5. As 
a consequence, the SATANS code has been modified to incorpor- 
ate a more accurate representation of the fluid and its inter- 
action with the motion of the fuel tank walls. The main 
elements of this modification are presented in Refs. 5 and 6. 

Both the original SATANS code and the modified code do not 
predict the fluid pressure due to the penetrating projectile. 
This pressure field can be provided by the Naval Weapons Cen- 
ter (NWC) Hydraulic Ram Program — Version One (HRP-Vl) , Ref. 7. 
The incident fluid pressure at the fluid-wall interface ob- 
tained from this code is used as input data for the modified 
SATANS code. For the very-early-time , short-duration shock- 
phase pressure, the Yurkovitch theory can be used. Ref. 8. 

An indepth hydraulic ram-structural response study per- 
formed by California Research and Technology, Incorporated 
(CRT) , Ref. 9, has been chosen for a comparison with the 
results obtained using the modified SATANS and HRP-Vl codes 
in order to develop confidence in both models. CRT compared 
their computer solutions with an experimental study performed 
by the University of Dayton Research Institute (UDRI) for veri- 
fication of their fluid and tank wall response predictions. 

The CRT study used a two-phase analysis with water as the fluid 
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and a thin circular plate for the entry wall. In phase I, 
the AFT0N-2D finite difference code was used to calculate the 
penetration dynamics of a sphere and to determine the fluid 
pressure at the tank wall. In phase II, the tank wall elastic- 
plastic response was calculated using the N0NSAP-2D finite ele- 
ment code. The pressure at the wall computed in phase I was used 
as the loading condition for the finite element analysis. Sec- 
tion II of this thesis presents the results from both studies 
for the projectile velocity and position in the fluid as a‘ 
function of time, the wall displacement as a function of time, 
and the circumferential stress vs. radius from center plate at 
specific times. A comparison of the crack damage prediction 
with the CRT results is made in section III using a crack damage 
prediction model developed at NPS by Fahrenkrog, Ref. 10. 
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II . COMPUTER STUDIES 



A . BACKGROUND 

When a fluid-backed plate is subjected to pressure from 
an external source or from within the fluid, the resulting 
dynamic response of the plate is signficantly affected by the 
contained fluid. The pressure incident upon the plate causes 
it to deflect and the fluid pressure at the fluid-plate inter- 
face is dependent upon the plate velocity, which leads to a 
coupled relationship between the fluid pressure and the plate 
motion. Ref. 11 contains an extensive review of this inter- 
action theory and of some of the finite element and finite 
difference computer codes used in the analysis of fluid-struc- 
ture response. 

One procedure that eliminates the coupled process is the 
plane wave piston theory. Fluid-structure interaction has 
been analyzed using this theory as early as 1940, Ref. 12. 
Piston theory acknowledges that the pressure radiated from the 
plate due to plate motion varies both in the transverse plane 
and in the direction normal to the plate; however, it ignores 
the transverse pressure variation in the equation governing 
the pressure field. Ball and Stronge, Ref. 5, have shown that 
piston theory applied to a water medium does not provide a 
sufficiently accurate solution for plate motion except at very 
early times. Their comparison of plate motion, except based 
upon piston theory, with both a NASTRAN solution and a modified 
SATANS solution revealed that the influence of pressure relief 
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in the transverse direction due to non-planar plate motion 
was significant and should not be neglected. 

As a result of the above study SATANS was modified by Ball 
to incorporate a more accurate representation of the fluid 
and its interaction with the containment vessel walls. A sub- 
routine was prepared that computes the fluid pressure through- 
out a cylindrical volume due to tank wall motion. The approach 
used was to approximate the equations of motion of the fluid 
by finite difference equations, Ref. 6 and Fig. 1. 

B. ENTRY WALL RESPONSE-A COMPARISON WITH CRT RESULTS 

To verify the modified SATANS code an indepth study per- 
formed by California Research and Technology was chosen for 
comparison with the SATANS/HRP-Vl results. The problem con- 
sidered by CRT is shown in Fig. 2. The CRT study utilized a 
two-phase analysis. In phase I, the penetration dynamics of 
an 11.1 mm diameter spherical penetrator and the time history 
of the fluid pressure at the fluid-structure interface were 
calculated from an adaptation of the AFT0N-2D finite differ- 
ence code. Pressure relief at the fluid-structure interface 
was accounted for by mass or inertia modeling of the entry 
wall. Thus, the fluid response is restricted only by the 
plate mass, and the effect of the plate elasticity is neglected 
in this phase. The fluid pressure computation in this phase 
was carried for 178 ysec. In phase II, the entry wall response 
was calculated using the NONSAP finite element code. NONSAP 
is a three-dimensional code that accounts for both material 
and geometrically nonlinear behavior. The interface fluid 
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pressure computed in phase I was used as the fluid loading on 
the entry wall. The wall response was computed for time much 
larger than the 178 ysec. for which the fluid pressure was 
calculated. For the larger time, the fluid pressure on the 
wall was neglected. Thus, the plate was essentially vibrating 
in a vacuum after 178 ysec. 

For the computer solutions using SATANS/HRP-Vl , the tank 
dimensions were assumed to be 200 in. x 200 in. x 100 in. deep, 
with the circular entry panel symmetrically centered on the 
tank front wall. The dimensions were selected so that during 
the period of the test runs there would be no reflected images 
from the other surfaces to interact within the fluid at the 
plate interface. Input data for HRP-Vl was assembled following 
the guidelines in the HRP-Vl user's manual. Ref. 7, with the 
code modified for use on the NPS IBM 360 system. 

Ballistic data for the penetrating sphere are presented 
in Fig. 3. The sphere was assumed to pass directly through 
an 11.1 mm diameter hole in the entry panel and enters the 
fluid at coordinates 100,100,0. The HRP-Vl code was used to 
obtain the entry wall pressures for the first 178 ysec. The 
fluid interface pressure was not computed beyond 178 ysec in 
order to make a comparison with the CRT results. Finite differ- 
ence mesh spacing for the SATANS/HRP-Vl analysis was 0.738 in., 
and the time step used was 2 ysec. 

The displacement and velocity of the spherical projectile 
calculated as a function of time by the AFTON-2D and HRP-Vl 
computer codes are presented in Figs. 4a and 4b, respectively. 
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An examination of these two figures reveals excellent agree- 
ment between the results from the two codes. According to 
Ref. 1 the pressure predicted by these two codes is also in 
good agreement with the experimental data obtained in the 
UDRI tests . 

Four different computer runs were made for the entry wall 
response using the SATANS/HRP-Vl codes. Two of the runs con- 
sidered a linear structural response; one run had the fluid- 
structure interface pressure relief due to wall motion accounted 
for, and the other run did not consider the pressure relief. 

The remaining two runs considered geometrically non-linear 
structural response, again, one with fluid-structure interface 
pressure relief accounted for, and one without. All four runs 
assumed no applied pressure beyond 178 ysec. Figure 5a contains 
entry displacement at 200 ysec. vs. distance from the impact 
point obtained from the four SATANS runs and from the CRT study. 
Figures 5b through 5j present the SATANS geometrically non- 
linear wall response, with interface pressure relief accounted 
for, at 200 ysec time intervals up to 2000 ysec. The corres- 
ponding data from the CRT study are also presented in these 
figures. 

Examination of Fig. 5a reveals that the effects of geo- 
metric non-linearity are small compared to the neglect of 
interface pressure relief, and that the SATANS/HRP-Vl results 
predict smaller displacement than does the CRT procedure. 
.However, note that the SATANS/HRP-Vl results do not include 
the effects of the initial shock pressure due to the initial 
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penetration predicted by the Yurkovitch theory. Ref. 8. In 
addition, SATANS results were computed for a plate that did 
not have an 11.1 mm hole at the plate center, and was linearly 
elastic. The neglect of the shock pressure probably leads to 
an underestimation of the entry wall displacement. The 
presence of the plate where the hole should be accounts for 
the zero slope at the plate center and probably further re- 
strains the plate motion, as does the linearly elastic 
restriction. 

Examination of Figs. 5b- j reveals that the maximum dis- 
placement beyond the hole predicted by the two codes is in 
very close agreement. Looking at the figures sequentially 
reveals disagreement, but still a good agreement in relative 
magnitude. 

Figures 6a-d present the circumferential stress vs. dis- 
tance from the plate center at 50, 100, 150, and 200 ysec for 
both codes. Examination of these figures reveals that the 
SATANS data are in excellent agreement with the CRT data in 
predicting the radial distance beyond which the circumferen- 
tial stress will be less than the nominal 70 kpsi tensile 
stress of the 2024-T3 aluminum plate. In addition, the rela- 
tive magnitudes of the predicted stresses are in good agree- 
ment beyond the point at which the stress drops below the 70 
kpsi level, however, there is a phasing disagreement similar 
to that described for Figs. 5b- j . Except for the stress at 
50 ysec, the maximum stress predicted by SATANS is higher than 
that predicted in the CRT study. 
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C. EXIT WALL RESPONSE-A COMPARISON WITH HRP-Vl MOD RESULTS 

The Naval Weapons Center, in the course of studies of the 
hydraulic ram phenomenon, developed a modification to their 
hydraulic ram pressure prediction code, HRP-Vl, that computes 
the response of a square exit plate to the pressure predicted 
by HRP-Vl. This modified code, known as HRP-Vl MOD, is ex- 
tremely quick in execution when compared to the considerably 
more complex codes like SATANS or BR-1HR. 

The original pressure code was modified to include compu- 
tation of the exit wall response by incorporation of the sub- 
routine VMIRROR. This subroutine provides a variable mirror- 
image response about the exit wall based on a pressure forcing 
function that allows incorporation of the pressure interac- 
tions into the description of the wall deflection. The spe- 
cific response model used replaces the free-surface mirror 
image for the exit wall used in HRP-Vl with an infinite void. 

A pressure distribution and a fluid velocity are calculated 
at the exit wall. The structural response of the exit wall 
is based on the assumption that the plate is simply supported 
and deflects in a double sine wave form. Both bending and 
stretching of the middle surface are accounted for. The inci- 
dent pressure distribution on the plate is calculated from 
the pressure due to the fluid-projectile interaction. From 
piston theory, a plane reflected wave is assumed which propa- 
gates parallel to the normal to the exit wall. Given initial 
pressures and fluid velocities at specific time increments, 
the required pressures and time varying parameters can be 
calculated and used to predict the exit wall deflections. 
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The problem selected for comparison with HRP-Vl MOD is 
shot 14 of the NWC test program, as reported in Ref. 8. The 
input data for both SATANS and HRP-Vl MOD are given in Ref. 8. 
The output data from SATANS and HRP-Vl MOD for the deflections 
at the center of the exit plate as a function of time are pre- 
sented in Fig. 7. For the time period up to 200 ysec, exami- 
nation of Fig. 7 reveals that the maximum center plate deflec- 
tions are in good agreement with a difference of approximately 
8%. A plot of the SATANS data for the plate deflection along 
a radius and the HRP-Vl MOD assumed deflection at 200 ysec is 
presented in Fig. 8. Relative magnitudes of the deflections 
compare favorably as does the phasing relationship. 

D. CRACK DAMAGE STUDY 

A simple procedure for estimating the length of a crack 
in a fuel tank wall due to hydraulic ram was presented by 
Fahrenkrog in Ref. 10. Fahrenkrog's procedure is based on 
several assumptions and observations. First, the fracture 
surfaces of tank wall plates tested appeared similar to the 
fracture surfaces on thin plate tensile test specimens, and 
therefore tensile test data from the stretching of cracked 
plates can be used. Second, there is a critical stress level 
for each crack length below which cracks in the tank walls 
will not propagate. 

The CRT case A study and UDRI tests provided data which 
can be used to substantiate Fahrenkrog's hypothesis about a 
critical stress level for crack propagation. The calculated 
circumferential stresses in the entry wall presented in Figs. 
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6a through 6d reached a maximum close to the plate axis at 
very early times (less than 200 ysec.) and then decayed 
rapidly with increasing radius. The stresses above the nomi- 
nal 70 kpsi tensile stress of 2024-T3 aluminum extend out to 
a radius of about 4 cm. 

The CRT case A corresponds to UDRI test FT5B conditions. 
Ref. 1. Examination of this test panel revealed that although 
the test panel contained inelastic distortion, radial cracks 
did not propagate outward from the entrance hole past about 
4 cm. 

Figure 9 is a graph which contains Fahrenkrog's plot of 
stress required for crack propagation versus cracklength from 
Ref. 10 with results for maximum circumferential stress versus 
radius from the SATANS analysis at t=200 usee. Note that the 
two curves intersect at a point slightly less than 5 cm., and 
at a stress level of nearly 39 kpsi. Briefly stated, this 
means that cracks that form around the entry hole should not 
propagate beyond this length at this time, which is in very 
good agreement with the UDRI test results. 

Fahrenkrog's maximum crack length prediction criterion 
requires the comparison of the maximum circumf erential stress 
along the radius of the plate with the stress required for 
crack propagation. A study of the stresses predicted by SATANS 
for times somewhat larger than 200 usee revealed that the 
stress levels dropped, indicating that the crack would not 
propagate beyond the length shown in Fig. 9. However, for 
much longer times, i.e, 1 msec — t — 2 msec, the stress levels 
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rebounded and peak above those shown in Fig. 9 for radii 
larger than 4 cm. These maximum stresses are shown in Fig. 
10. Hence, a strict application of Fahrenkrog's criterion 
in Fig. 10 leads to a predicted maximum crack length of 18 
cm. 
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III. CONCLUSIONS AND RECOMMENDATIONS 



A. ENTRY WALL RESPONSE 

The NWC HRP-Vl computer code is an excellent tool for pre- 
dicting the fluid pressure field caused by a projectile pene- 
trating a rectangular prism of fluid. When coupled with the 
current modified SATANS code, the two codes provide a good means 
of predicting the response of the entry wall of a fluid-filled 
tank subjected to a penetrating projectile. However, during 
this study the SATANS code was unable to accurately predict 
the response of an entry panel when a central hole was intro- 
duced to allow projectile entry directly into the fluid. 

The modification of the subroutine OUTPUT in the existing 
SATANS code is an invaluable tool for the analyst because it 
provides the test plate stresses or strains at discrete radial 
locations at various times . 

B. EXIT WALL RESPONSE 

SATANS and HRP-Vl MOD provide close results when predicting 
the exit wall response at the plate center for early times. 
Furthermore, the HRP-Vl MOD code uses approximately 1/3 the 
computer time required by SATANS to predict the response of the 
exit plate under the same conditions. However, the HRP-Vl MOD 
code assumes a sinusoidal plate deflection which, when compared 
to deflections predicted by SATANS, gives a plate displacement 
that may be too smooth. 

At the present time, the HRP-Vl MOD program is limited to 
predicting only exit wall response. Further study should be 
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performed with the HRP-Vl MOD code to determine if it can be 
modified to also predict entry wall response. 

C. CRACK DAMAGE STUDY 

The SATANS code with the modified OUTPUT subroutine prints 
stress levels in the entry panel which are in reasonable 
agreement with the CRT and UDRI test results for stresses be- 
low the maximum tensile stress, up to 200 ysec. When the 
SATANS-predicted circumferential stresses at 200 ysec, versus 
radius, is applied to Fahrenkrog's crack length prediction 
criterion, the predicted crack lengths are nearly the same as 
the lengths actually measured under controlled testing by 
UDRI test FT5B. When maximum predicted circumferential stress 
vs. radius is applied to Fahrenkrog's criterion, the maximum 
predicted crack length is appreciably longer. 

Because the SATANS data in conjunction with Fahrenkrog's 
crack length prediction criterion predicted conservative 
results compared with the actual test results from UDRI, fur- 
ther comparison should be performed with actual test results 
from several different tests. These comparisons will help to 
establish a standard for crack length prediction ij. fluid- 
filled tanks subjected to a penetrating projectile. 
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